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The work described in this report was performed under 
Contract No. NAS2-10765 for the National Aeronautics and Space 
Administration, Ames Research Center. The technical monitor for 
NASA-Ames was Dr. William Warmbrodt. The principle author from 
Systems Control Technology, Inc. was Mr. Russell L. Hull. Mr. 
James H. Vincent was program manager. Technical support was 
provided by Mr. James W. Fuller. Programming support was 
provided by Mr. Bob Bullock. Report preparation efforts were 
directed by Ms. Clare Walker. 
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NOTATION 


ail, ai2. 

321, 

322 = components of the A dynamics matrix 

A 

= 

* 

coefficient matrix in state equation x = Ax + Bu 

AEE 

= 

engine state coefficient matrix in engine state 
equation 

Arfu 

= 

fuselage state coefficient matrix in rotor state 
equation 

Arr 

= 

rotor state coefficient matrix in rotor state 
equation 

Art 

s 

drive train state coefficient matrix in rotor 
state equation 

ATEl, ATE2 

= 

engine 1 and 2 state coefficient matrix in drive 
train state equation 

ATT 

= 

drive train state coefficient matrix in drive 
train state equation 

AlS 

= 

main rotor lateral cyclic control input 

A('i') 

= 

coefficient matrix in rotating frame blade 
flapping equation 

A'(4;) 


coefficient matrix in rotating frame blade 
lagging equation 

t>l, b2 

- 

components of the B input matrix 

B 

3 

coefficient matrix in state equation x « Ax + bu 

BEF 

= 

fuel control output coefficient matrix in engine 
state equation 

BET 

3 

drive train output coefficient matrix in engine 
state equation 

BETA 

= 

collective setting 

Brc 

3 

rotor control coefficient matrix in rotor state 
equation 

Brfu 

3 

fuselage output coefficient matrix in rotor state 
equation 
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NOTATION (Continued) 


Brl 

Brt 

BTL 

BTR 

Bis 

B(^) 

B'(«K) 

c/Cc 

Cl. C2, C3 

C 

Cafu 

Car 

CeE 

CpG 

CgP 

Crfu 

Crr 


= lag damper coefficient matrix in rotor state 
equation 

* drive train output coefficient matrix in rotor 
state equation 

= lag damper output coefficient matrix in drive 
train state equation 

= rotor output coefficient matrix in drive train 
state equation 

= main rotor longitudinal cyclic control input 

= coefficient matrix in rotating frame blade 
flapping equation 

= coefficient matrix in rotating frame blade 
lagging equation 

* damping ratio 

= coefficients 

= coefficient matrix in output equation y = Cx + Du 

= fuselage state coefficient matrix in fuselage 
aero output equation 

= rotor state coefficient matrix in fuselage aero 
output equation 

= engine state coefficient matrix in engine output 
equation 

= parameter, fuel flow to gas generator dynamics 

= parameter, gas generator to power turbine dynamics 

= fuselage state coefficient matrix in rotor output 
equation 

= rotor state coefficient matrix in rotor output 
equation 
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NOTATION (Continued) 


CrT 


drive train state coefficient matrix in rotor 
output equation 

CTEl, CTE2 

s 

engine state coefficient matrix in drive train 
output equation 

Ctt 

s 

drive train state coefficient matrix in drive 
train output equation 

Cl 

3 

speed ratio, tail rotor/main rotor 

de 

= 

equivalent viscous damping between engine and 
fuselage 

dH 

3 

equivalent viscous damping between rotor hub and 
fuselage 

dHE 

3 

equivalent viscous damping due to shaft rate of 
deflection between engine and hub 

D 

3 

coefficient matrix in output equation y * Cx + Du 

^AC 

= 

rotor control coefficient matrix in fuselage aero 
output equation 

Dafu 

= 

fuselage output coefficient matrix in fuselage 
aero output equation 

Dari 

3 

rotor accelerations coefficient matrix in 
fuselage aero output equation 

Dat 

3 

drive train output coefficient matrix in fuselage 
aero output equation 

DeF 

3 

fuel control output coefficient matrix in engine 
output equation 

DrC 

3 

rotor control coefficient matrix in rotor output 
equation 

Drfu 

3 

fuselage output coefficient matrix in rotor 
output equation 

Drri 

3 

rotor acceleration output coefficient matrix in 
rotor output equation 
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NOTATION (Continued) 


Drt 


= drive train output coefficient matrix in rotor 
output equation 


ECU 


= electronic control unit 


^o» flc> fl<?> f2c» f2«; = components of matrix 

expaSled in a truncated series 

Pi = i^h force component 


F(^) 


coefficient matrix in rotating frame blade 
flapping equation 


F'(^) 


coefficient matrix in rotating frame blade 
lagging equation 


G( 


coefficient matrix in rotating frame blade 
flapping equation 


G'(^) 


coefficient matrix in rotating frame blade 
lagging equation 


H 


rotor forward force 


H 

HMU 

'g 

Ih 

J 


* transpose of transformation matrix from rotating 
blade coordinates to multi-blade coordinates 

= hydroraechanical unit 

= inertia of gas generator turbine and compressor 
= inertia of rotor hub 

= inertia of power turbine and output shaft 
= matrix 


J 

^Hl, kH2 


IDS 

Lp 

Lh 


= rotor lateral force 

= equivalent torsional stiffness between the rotor 
hub and engine 1 and 2 

* load demand spindle 

= fuselage roll aerodynamic torque 

= rotor roll torque 


X 



NOTATION (Continued) 


Mp 

Mh 

NG 

NP 

Np 

"g 

"p 

Npi, Np2 
P 

PAS 

PO 

Pr 

PS3 

P4I 

P45 

q 

Qcomp 

Qe 

Qgturb 

Qh 

Qld 

Qmr 


= mass 

= fuselage pitch aerodynamic torque 
= rotor pitch torque 
= gas generator turbine speed 
= power turbine speed 
= fuselage yaw aerodynamic torque 
= gas generator turbine speed (rpm) 

= power turbine speed (rpm) 

* shaft speed of power turbines 1 and 2 

* roll rate, fuselage 

» power available spindle 
= ambient pressure (altitude) 

= compression ratio 

a compressor output static pressure (psi) 
a gas generator turbine inlet pressure (psi) 
a power turbine inlet pressure (psi) 
a pitch rate, fuselage 
a compressor load torque 
a engine output torque (ft lb) 
a gas generator turbine output torque 
a rotor yaw torque 
a sum of lag damper torques 
a main rotor torque 
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NOTATION (Continued) 


Qpturb 
Qpi. QP2 

Qrotor 
Qi» Qz 

r 

^Hl, th2 


Rl 

s 

SHP 

T 

Tin 

Tr 

Tz 

T45 

u 

u 

V 

VG 

w 

W4SR 

Wp 


= power turbine output torque 

* torque transferred by drive train to power 
turbines 1 and 2 (ft lb) 

* helicopter load torque on engine 
= torque, engines 1 and 2 

* yaw rate, fuselage 

= gear ratio between the main rotor hub and engines 
1 and 2 (rotor speed/engine speed) 

= variable rate limit 

a Laplacian operator 

= shaft horse power 

= final time 

= rotor axial force 

= temperature in 

* engine thrust before rate limiting 
> temperature, ambient 

» power turbine inlet temperature 
= axial velocity, fuselage 
= input 

* lateral velocity, fuselage 
= variable guide vane setting 
= vertical velocity, fuselage 
= engine mass flow condition 
= fuel flow rate (pps) 


xii 



NOTATION (Continued) 


X 

= 

vector 

of 

states 


= 

vector 

of 

engine states 

XI, X2 

3 

states 

1 

and 2 

^El. Xe2 

= 

vector 

of 

states for engine 1 and engine 2 

Xp 

3 

fuselage 

forward aerodynamic force 

Xt 

= 

vector 

of 

drive train states 

y 

3 

vector 

of 

outputs 

XE 

= 

vector 

of 

engine outputs 

yp 

3 

vector 

of 

fuel control outputs 

yi 

3 

vector 

of 

lag damper outputs 

yp 

3 

vector 

of 

rotor outputs 

yT 

3 

vector* 

of 

drive train outputs 

Yp 

» 

fuselage 

lateral aerodynamic force 

Zp 

a 

fuselage 

downward aerodynamic force 

GREEK, ETC. 





a 

= 

linear 

parameter 

6 

s 

vector 

of 

multi-blade flapping coordinates 

Sr 

a 

vector 

of 

rotating frame flapping coordinat 


®o» ^lc» ®ls» ®2c» ®2s “ multi-blade flappine 

coordinates for a five-bladed rotor 

®1> ®2» ®3» ®4» ® rotating frame flapping 

coordinates tor a five-oladed rotor 

y = ratio of specific heats 
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NOTATION (Continued) 


*^7 = power lever angle 

^ 1 , ^2 ® deflections in the drive train between engine 1 

and the rotor hub, and engine 2 and the rotor hub 
(referenced to the hub frame) 

C * vector of multi-blade lagging coordinates 

Cr a vector of rotating frame lagging coordinates 

^o» ^lc> ^ls» ^2c» ^2s “ multi-blade flapping 

coordinates for a five-bladea rotor 

^1» ^2» ^3» ^4» ^5 ■ rotating frame flapping 

coordinates for a five-oladed rotor 

* compressor efficiency 

® * roll attitude, fuselage 

®c = main rotor collective control input 

®TR = tail rotor collective control input 

u * rotor advance ratio 

I = summation signal 

= time constant 


"P 

% 

^MR 


* gas generator turbine dynamic time constant 
= power turbine dynamics time constant 
= roll attitude, fuselage 
= yaw attitude, fuselage 
= azimuth angle of i^h blade 
= power turbine rotation rate 
= main rotor rotation rate 
= tail rotor rotation rate 


xiv 



NOTATION (Concluded) 


SUBSCRIPTS, SUPERSCRIPTS, ETC. 


( he 

( )n 
( )o 

( )o 

( )trim 

A( ) 

5 ( ) 

*( ) 

/( )/ 

( )T 
( )-l 
(■) 

(') 


( ) at initial condition 

new ( ), refers to condition at current time 
old ( ), refers to condition at previous tim< 
trim or initial value of ( ) 

trim condition of ( ) 

perturbation of ( ) from trim 

perturbation of ( ) 

the axis transformation components have been 
deleted from ( ) 

absolute value of { ) 

transpose of ( ) 

inverse of ( ) 

vector ( ) 


step 

step 
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I. INTRODUCTION AND SUMMARY 


1.1 INTRODUCTION 

It is important to be able to predict coupled propulsion 
system and rotor system dynamic responses for rotorcraft. The 
capability is needed for assessing the effects of varied fuel 
control design, rotor lag damper characteristics, and drive 
system equivalent flexibility, among other possible 
improvements. Specific occurrences of coupled propulsion 
system/rotor system dynamic problems have been observed in 
current rotorcraft, such as the CH-47C. The potential exists for 
stability problems in new vehicles, such as the UTTAS, AAH, and 
XV-IS, and in prototype vehicles designed for enhanced 
performance, such as the X-Wing or HH-20 (CCR). Experimental and 
operational data from rotorcraft indicate that improper matching 
of rotor dynamic characteristics with drive system and engine 
response modes is a fundamental source of undesirable vehicle 
stability and vibration characteristics. The technology applies 
to a wide range of rotor system types and propulsion/fuel control 
configurations . 

Development of the capability to predict coupled 
propulsion/rotor dynamic responses for a range of rotorcraft 
designs has not been systematically initiated. The basic reason 
is that the development of a comprehensive analysis approach 
using a generic methodology, is a highly multidisciplinary 
objective. The objective is broader than the specific product 
development and manufacturing objectives of many helicopter or 
engine companies. 
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1 . 2 SUMMARY 


Systems Control Technology, Inc. (SCT) undertook a program 
to establish a generic methodology and model structure for 
performing coupled propulsion/rotor response analysis that would 
be applicable to a variety of rotorcraft types. In this program, 
SCT performed the following tasks: 

(1) Developed a model structure, adaptable to a wide range 
of rotorcraft configurations, for simulating coupled rotor/ 
propulsion dynamics. 

(2) Defined a methodology for parameterizing the model 
structure to represent a particular rotorcraft. 

(3) Constructed a nonlinear coupled rotor/propulsion model 
for a particular (though not existing) rotorcraft as a test case 
to use in analyzing coupled system dynamics. The model included 
a propulsion system and a rotor/fuselage system as follows: 

(a) The propulsion system was constructed from an 
existing engine/fuel control model and a drive train 
model derived by the author. 

(b) An existing helicopter rotor and airframe handling 
qualities simulation was used for the rotor fuselage 
system . 

(4) Constructed a partially linear coupled rotor/propulsion 
model by applying perturbation and subset regression techniques 
to linearize the engine and rotor components of the nonlinear 
model. This step was not completed due to difficulties 
associated with the periodic dynamic effects in the rotor 
dynamics at the rotorcraft speed chosen for linearization. 

1.3 ORGANIZATION OF REPORT 

This report is divided into two volumes. Volume I contains 
the body of the report and the analytical part of the appendix. 
Volume II contains documentation on the computer models that were 
developed . 
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In Volume I, Chapter II gives a problem overview and 
discusses the background of the problem. Chapter III describes 
the generic modeling methodology. Chapters IV and V document the 
development of the propulsion system and the rotor/fuselage 
models, respectively. Chapter VI describes the formulation of 
the resulting coupled rotor/propulsion system model and describes 
a test case that was developed. Finally, Chapter VII summarizes 
the project and makes recommendations for further model 
development . 
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II. BACKGROUND AND PROBLEM OVERVIEW 


This section reviews the background of the program and 
details the issues associated with development of a coupled rotor 
and propulsion system analysis. 

2.1 BACKGROUND AND PROBLEM OVERVIEW 

2.1.1 Background 

r 

The basic dynamic elements of a typical helicopter 
propulsion drive system are shown in Figure 2.1. One or more 
engines drive a system of shafts, gear reductions, the main 
rotor(s) and the tail rotor. The main and tail rotor are the 
primary force and moment producers. While many elements of this 
sytem have been designed to minimize coupling between the 
elements, recent experience has demonstrated the need for a more 
systematic treatment of rotor/propulsion system interaction. 

In older helicopters, the flight and propulsion system 
dynamics of the rotorcraft were not closely coupled because there 
was a separation between the speeds of response of these 
subsystem dynamics. Propulsion system dynamics were much faster 
than the flight dynamics, and thus handling qualities designers 
could assume the rotor had a constant rate as far as handling 
qualities were concerned. In recent helicopter designs, however, 
demands for higher flight performance have led to designs in 
which there is significant and even destabilizing coupling 
between subsystems. 

One dramatic example of this involved the CH-47C helicopter 
with a Lycoming T-53 engine [1]. During testing of this vehicle, 
severe rotor oscillations were encountered, leading to pronounced 
fuselage vibration. Time histories of the responses of several 
of the variables that were involved are shown in Figure 2.2a 
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Figure 2.1 Typical Helicopter Rotor and Transmission System 
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Figure 2.2a Torque Oscillation Flight Test Data [1] 



Figure 2.2b Lag Damper Force-Velocity Curves [1] 
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Figure 2. 2d Final Torque Oscillation Simulation [1] 
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[1]. Two parameters were varied in the test program to reduce 
the oscillation shown. First, the blade lag damper was "tuned 
down" which reduced that particular problem (Figure 2.2b), and 
then the fuel control time constant was increased and the gain 
reduced (Figure 2.2c). The result was the acceptable response of 
Figure 2. 2d. 

In response to the CH-47C problem, and related experiences, 
a series of studies was sponsored by the Army at Ft. Eustis. 
Results of these studies are detailed in Ref. 2. Table 2.1 lists 
some of the relevant interface problems that are discussed in the 
volumes describing these studies. 

The background for the current program was the recognized 
need for a systematic approach to analyze and avoid such problems 
in design and test programs. 

2.1.2 Problem Overview 


The principal elements which must be considered in a 
systematic rotor/propulsion system response analysis are: 

• engine models; 

• fuel control models; 

• rotor model; and 

• drive train models. 

The frequency range of interest initially is a range of 
approximately 0-5 Hz (or zero to approximately once per rotor 
revolution). A model valid over a larger frequency range would 
be desirable. However, this would require more detailed 
representations of the rotor system, drive train, and aircraft. 
This effort is beyond the scope of the present analysis and 
consequently application of the model is restricted to this 
frequency range. Still, a significant number of previous 
engine/airf rame/dri ve train dynamic interface development 
problems [2] have occurred within this frequency range, including 
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Table 2.1 


Interface Problem History 


• TORSIONAL MODE OSCILLATION - COLLECTIVE LAG 

SIKORSKY CH-53E/T64 

BELL 0H-4A/D0A 250 SERIES 

BOEING CH-47C/LYC0MING T55-L-n 

• EXCESSIVE ROTOR DROOP i OTHER RPM CONTROL PROBLEMS 

SIKORSKY/BLACKHAWK 
BELL YAH-63/GE T700 
BOEING YUH-61A/GE T700 

• ENGINE HUNTING 

SIKORSKY XH-59A/PT6T-6 

• ENGINE LOAD SHARING OSCILLATION - 2 ENGINES 

SIKORSKY SH-34H/GE T5 
TILT ROTOR 

• EXCESSIVE TRANSIENT LOADING (HIGHER FREQUENCY) 

BELL YAH-63/GE UOO - OVERSPEED TRIP 
BOEING CH-47/LYC0MING T55 - ROTOR START UP 
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engine torque oscillations, engine load-sharing problems, and 
excessive transient RPM droop. 

Modern rotorcraft utilize free-turbine engines where the 
second or power turbine section is connected to the output shaft 
(Figure 2.3 [3]). Because the load on this output shaft varies 
with collective control and flight speed, modern turbines, in 
particular, require multiple sensing feedback paths and 
multivariable controls to sustain design RPM. The control system 
of the current GE/T700 engine, for example, is illustrated in 
Figure 2.4 [3]. 

The propulsion system dynamics are especially important to 
flight control in the vertical mode because the propulsion system 
is then the only source of power. To change the amount of lift 
over even a few seconds, the fuel flow must be changed. In 
forward flight, kinetic energy can be exchanged for lift. In the 
vertical mode there is no mechanism of energy storage that is 
large enough to provide a sustained engine output variation. 
Because of the close relationship between flight and engine 
performance, overall performance can be enhanced by jointly 
controlling these modes. 

In current practice, the engine and flight control system 
are designed separately. Handling quality analyses and automatic 
flight control design analyses do not typically consider the main 
rotor RPM a degree of freedom. In addition, specific goals for 
engine transient response are not well defined. Instead, the 
engine is specified to respond to a step input as quickly as 
possible without overly degrading stability margins. This 
approach has its limits, especially when the engine/fuel 
controller response is not fast enough for good handling, but is 
so fast that it interacts with some torsional/f lexural modes and 
causes vibrational oscillations. This situation can often occur 
in rotorcraft design because of the relatively low-frequency 
torsional bending mode of the rotor. Figure 2.5 [4] shows some 
of the torsional modes for the XV-15. The first mode is shown to 
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0. FREE STREAM AIR CONDITIONS 

1. INDUCTION SYSTEM OUTLET/ENGINE INLET INTERFACE 

2. COMPRESSOR INLET 

3. COMPRESSOR OUTLET 

3.1 BURNER INLET 

4. BURNER OUTLET/FIRST TURBINE STAGE INLET INTERFACE 

4.1 FIRST TURBINE ROTOR STAGE INLET 

4.5 GAS GENERATOR TURBINE OUTLET/POWER TURBINE INLET 
INTERFACE (INLET TO POWER TURBINE NOZZLE) 

5. LAST TURBINE STAGE OUTLET 

7. ENGINE OUTLET/EXHAUST DUCT INLET INTERFACE 
9. EXHAUST DUCT OUTLET 


Figure 2.3 Elements of Current Rotorcraft Engines [ 
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Figure 2.4 T700 Control System [3] 
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Torsional modas 


Fraquency (per rev) 


Mode shapM 


1st asymmetric flexible mode 
(interconnect shaft mode) 

56S rpm 

458 rpm 

Pylor 

Roto 

r Interconnect See Note 

39 

48 

1st symmetric flexible mode 
(symmetric power turbine mode) 

1 26 

1 S5 

' 

n 


2nd asymmemc flexible mode 
(asymmetric power turbine mode) 

1 43 

1 79 

/ 



2nd symmetric flexible mode 
(symmetric pylon roll mode) 

3 42 

4 22 

1 

1 

w 

3rd asymmetric mode 
(asymmetric pylon roll) 

3 49 

430 

y 




Note Engine rotations referred to tilt rotor 

Positive engine rotation same as corresponding tilt rotor 


Figure 2.5 


Torsional Rotor/Drive Train/Engine Bending 
Modes for XV-15 [4] 
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be at 3.67 Hz. In addition, many conventional articulated 
helicopters need to incorporate a lag damper to prevent ground 
resonance. This adds another set of torsional modes, the first 
of which involve rigid blade collective lag about the lag hinge, 
at approximately 3 Hz. The demand for more responsive 
helicopters in recent years has been met by using more responsive 
engines. However, due to torsional mode interaction, the 
limitation of this approach has been reached, and manufacturers 
have started to use combined automatic flight/propulsion control 
systems . 

In conventional helicopters, the propulsion system dynamics 
can cause handling problems and reduced responsiveness. A common 
mode of flight/propulsion interaction is known as rotor droop. 
Conventional rotorcraft regulate the rotor RPM to be near a 
design value which is limited in part by consideration of various 
rotor blade bending modes. When there is an external disturbance 
or a step-collective input, the increased drag on the rotor makes 
it slow down or droop. The increased torque demand is apparent 
to a separately designed fuel controller only through RPM droop, 
and therefore the engine response is delayed. This droop is 
caused by the transfer of rotor kinetic energy into lift. 

However, as the RPM decreases so does the lift and the 
effectiveness of the collective input. This degrades the 
handling quality of the vehicle. Also, when the rotor torque 
increases due to fuel controller action, the yaw reaction torque 
on the fuselage increases. This causes a perturbation in yaw 
that is further compounded when the main and tail rotor RPMs 
droop causing a reduction in tail rotor effectiveness. Pilots 
often compensate for droop by applying some pedal input along 
with collective stick input. The cross-coupling between the 
propulsion system and the yaw mode can also act in reverse. A 
change in yaw moment causes a change in the reaction torque on 
the rotor, which will cause an RPM variation. A fast fuel 
controller can detect this and cause the engine to change its 
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torque output. However, due to dynamic phase lag, this torque 
change acts not so much to regulate RPM, but rather to reduce the 
stability in the dutch roll mode. These fuselage/engine 
cross-coupling phenomena are discussed for the Blackhawk 
helicopter in Ref. 5. 

The XV-15 has some additional engine/airframe interactions 
that have caused trouble [6], The aerodynamic damping of the RPM 
governing system decreases with decreasing power required. This 
has caused the governor to oscillate in high-speed flight. The 
current solution is to lower the governing gain for good 
operation with worst-case damping. 

The state of the art solution to engine/airframe interaction 
is to keep the fuel flow controller slow enough to eliminate 
unfavorable coupling and to use a rudimentary form of combined 
control. The combined control involves using the collective 
stick input to bias the RPM error signal. In this way the engine 
is made aware of the changed torque requirement when the pilot 
commands it. This was used in the Blackhawk and, to a much 
greater degree, in the tilt rotor. When the rotor inertia is 
small, the link between the collective input and throttle control 
needs to be more pronounced. In the XV-15 tilt rotor, for 
example, the lever in what is normally the collective stick 
position is the power lever which is the engine throttle 
control. The collective rotor blade pitch control levers are 
cam-linked to the power lever. Also, the RPM is not regulated 
using fuel flow, but rather by adjusting the collective pitch of 
the rotors. This method of regulation prevents the exchange of 
rotor kinetic energy for lift. It is necessary because the 
rotors have relatively low moments of inertia and therefore do 
not store much energy. Exchange of energy causes larger RPM 
variations in smaller rotors. 

The helicopter or XV-15 response could be further improved 
if a controller were designed based on the combined coupled 
airf rame/propulsion system characteristics. 
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Another important engine control problem occurs with 
multiple engines. This problem is load-sharing oscillations, as 
demonstrated in Figure 2.6 [7]. An initial load-sharing problem 
(now solved) occurred in the Sikorsky SH-34H/T-58 helicopter. 

The major tool required for combined propulsion/flight 
control is a suitable dynamic model. 

2.2 METHODS OF ANALYSIS 

In the past, the procedure for designing for adequate 
combined system (engine, rotor, fuel control, etc.) performance 


has been 

to : 


(1) 

design the engine control for optimal response to 
step-load change; 

(2) 

design the rotor system assuming it 
by the engine at constant RPM; 

will be driven 

(3) 

use currently known rules of thumb 
the design for interactions between 
elements (e.g. rotor, fuel control. 

to account in 
system 
engine) . 


More recently, at Sikorsky for example [5], high-fidelity 
nonlinear models have been used to analyze the interactions of 
rotorcraft system components and to design to prevent or to 
correct interaction problems. Other approaches have involved ad 
hoc modifications to nonlinear models to allow consideration of 
specific interaction problems. 
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III. GENERIC MATH MODELING 


The purpose of this project was to develop a generic 
methodology and a model structure for performing coupled 
rotor/propulsion response analysis. The methodology was desired 
to allow accurate consideration of the interactions of rotorcraft 
system components, in a specific frequency range, without the 
complexity of high-fidelity nonlinear models. This chapter 
describes the method of approach that was developed. The 
sections of this chapter that follow identify the requirements on 
the model structure and the approach developed to meet these 
requirements . 

3.1 PURPOSE 

A rotorcraf t/propulsion dynamics model would be used in 
applications like the following: 

• flight/rotor/propulsion control system research and 
development 

*• trade studies for defining specifiction requirements, 
e.g. ensuring acceptable coupling of rotor and 
propulsion dynamics 

• predicting flight characteristics for flight test 
planning and flight safety analysis 

• support of development flight testing in areas such as 
troubleshooting, control system revision, and data 
analysis 

These applications require an accurate model, but do not require 
the high complexity of some nonlinear engine and rotor component 
models that are available, but not coupled together. 
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3.2 MODEL STRUCTURE REQUIREMENTS AND SCOPE 


The intended uses of the model require that it accurately 
predict rotorcraft handling qualities, instabilities, 
oscillations, and transient performance within a specified 
frequency range. The intended uses also require that the model 
structure and variables relate directly to physical explanations 
so the results of an application of the model can be easily 
interpreted in terms of the structure. 

Certain effects should be considered in the model. These 
include torsional mode stability, engine load-sharing problems, 
RPM droop and effects on handling quality. 

Only a limited frequency range is required for the model. 
Figure 3.1, compiled by Sikorsky Aircraft [5], shows frequencies, 
representative of a medium-sized helicopter, from which it is 
apparent that the dynamic modes of interest in this project occur 
at frequencies between zero and the shaft speed, This 

frequency range includes the first torsional modes of the 
rotor/drive system and rotor/airframe, the drive system/fuel 
control "hunting" frequency, fuselage instabilities, and the 
rigid flap and lag motions of the rotor blades. 

3.3 GENERIC MODELING APPROACH 

The sections that follow give a general picture of a generic 
modeling approach that was developed that can be used to 
construct a coupled rotor/propulsion dynamic response model for 
investigating performance in the zero to once per rotor 
revolution frequency range. 

3.3.1 Model Structure 


The generic modeling approach entails construction of a 
total system model from linear or nonlinear state models of each 
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Figure 3.1 Modal Frequencies of Interest In Engine/Fuel 

Control Design [5] 
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subsystem component (e.g. fuselage, rotor, engine). This 
approach provides a flexible, adaptive model structure that can 
be used to simulate a variety of rotorcraft designs. The 
structure allows subsystem models to be added, deleted, replaced, 
or modified to meet a particular application. 

Figure 3.2 shows the generic structure chosen for the 
rotor/propulsion dynamics model. The figure shown is for the 
dynamics of a twin engine, single main rotor helicopter in hover 
or forward flight. This figure contains a block for each 
subsystem model. Connecting lines show the paths of interaction 
between subsystems. The subsystems are linked along these paths 
by physical principles. 

The structure used for each component in the coupled system 
rotorcraft model is shown in Figure 3.3. This input/output form 
is general, and allows independent subsystem components to be 
easily connected together into a system. 

The subsystem model components are derived from detailed 
nonlinear simulations of each component. These detailed 
simulations can be obtained from the component manufacturers 
(e.g. for engines and fuel controls) or currently available 
simulations (e.g. the GENHEL rotor, airframe, control system 
simulation). 

The control system components of the model can be allowed to 
contain important nonlinearities. The control systems are the 
lag damper, fuel control, and flight control. Nonlinear effects 
are important in the control models because of their effects on 
the dynamic response. Also the model may be used in the design 
and analysis of control systems and thus may require 
high-fidelity nonlinear control models. 

3.3.2 Linear Model Formulation Methodology 

Models for the linear components of the combined system 
model are formulated using perturbation techniques. With these 
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Figure 3.3 Input/Output Form Used For Modeling Subsystems 
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techniques, equations are constructed using perturbation values 
about trim (or nominal) operating conditions of the rotorcraft. 
The trim values can be added back onto the predicted response 
perturbations to yield the "total" response values. 

Several techniques are available for extracting these 
linear, simplified models from complex nonlinear simulations. 
These include: 

• state and derivative perturbation techniques; 

• subset regression techniques; 

• system identification techniques with time 
histories generated from a nonlinear simulation 
model . 

Additionally, equations can be derived or extracted analytically. 
3.3.3 Subsystems 


The subsystems of the generic model shown in Figure 3.2 
include the following: 

• rotor 

• fuselage 

• engines (2) 

• fuel controls (2) 

• drive train 

• lag damper 

• tail rotor 

Some of the details and assumptions which were used in developing 
the generic form for modeling each of these subsystems are 
described below. 
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Fuselage 

The fuselage was assumed to be rigid. It was modeled by 
nine states. These included axial, lateral, and vertical 
velocities (u, v, w) ; roll, pitch, and yaw rate (p, q, r); and 
roll, pitch and yaw attitude (<P, a, . 

Rotor 


The rotor was modeled using multi-blade coordinates. 

Any of three approaches could have been used. These 
included rotating frame coordinates, multi-blade coordinates, and 
tip-path-plane coordinates. 

In rotating frame coordinates, each blade is modeled 
separately. This representation makes intersystem (e.g. between 
rotor and fuselage) dynamic coupling analysis more complex and 
results in a loss in engineering insight. A fixed-frame 
representation is more amenable to intersystem dynamic coupling 
analysis. A fixed-frame representation, such as in multi-blade 
or tip-path-plane coordinates, using constant coefficient values 
is adequate for many rotor/drive train/engine coupling modes. 

Because the fixed frame representation offers better insight 
into what is happening and is less complex, it was used in the 
generic model structure. A multi-blade rather than 
tip-path-plane representation was used with the number of main 
rotor flapping states and lagging states each to be equal to the 
number of blades. This number of states could later be reduced 
if it were found to represent higher frequencies than required in 
the model. 

Lag Damper 

Helicopter lag dampers are very nonlinear, as shown by the 
typical response plotted in Figure 3.4 [7]. Because the lag 
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damper has an important effect on the rotor dynamics, it should 
be modeled as a separate, nonlinear element in the model. This 
approach was taken. It allows straightforward analysis of the 
sensitivity of rotorcraft performance to variations in the lag 
damper design. 

Drive Train 

Drive trains are normally modeled as systems of torsional 
springs, dampers, and inertias. Gearing ratios add some 
complexity to this, but standard models [1,5,12,13] are still 
linear. Nonlinear backlash effects were not considered. 


Engine and Fuel Control 


The engine and engine control system models were required to 
describe the dynamics and coupling to the rotor and airframe up 
to the rotor RPM frequency (3 to 5 Hz). A perturbational model 
was used for the engine model while a nonlinear simulation was 
used for the fuel control. The perturbational engine model was 
extracted from a nonlinear engine simulation. 


Tail Rotor 


Tail rotor dynamics could either be modeled explicitly, as a 
separate part of the model, or as a control acting directly on 
the dynamics of the fuselage, without separate tail rotor 
dynamics. Because this was the first construction of the model 
the simpler approach was taken. At a later date, the structure 
could be extended to allow explicit modeling of the tail rotor. 
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Control Systems 


A flight control system model was necessary in order to test 
"fly" the model. A control system in an existing helicopter 
simulation was used. 

The pilot controls included in the model were: 

• main rotor - collective pitch, longitudinal and 
lateral cyclic pitch 

• tail rotor - collective pitch 
3.4 MODEL ORGANIZATION 

The generic model structure was constructed in a modular 
fashion with modules as were shown in Figure 3.2. In this 
manner, the modules could be developed independently. This 
formulation gave the generic model structure flexibility because 
any module could be independently modified or replaced by a new 
module or one parameterized for a different helicopter. The 
detailed development of each of the modules of the model that was 
developed to test and demonstrate the structure are discussed in 
the chapters that follow. 
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IV. PROPULSION SYSTEM MODELING 


A propulsion system model of a particular engine, fuel 
control and drive train combination was constructed n the format 
of the generic model structure. The model was constructed n 
order to test and demonstrate the methodology being developed. 
The development of this propulsion system model is described in 
this chapter. 

The elements which make up the propulsion system are: 

• engine models 

• fuel control models 

• drive train models. 

This chapter describes the development of this propulsion 

system model. The results of the chapter are a fully 
parameterized computer simulation modeling a specific helicopter 
engine/fuel control/drive train combination. 

4.1 DEVELOPMENT OF ENGINE AND FUEL CONTROL MODELS 

4.1.1 Problem Overview 


The propulsion system can be described in mathematical terms 
in several ways. The choice of model form has an important 
impact on the usefulness of the resulting model. A consideration 
in developing the model is the interface between the engine and 
the input/output boundaries. Models that describe the throttle 
linkage and fuel control as well as the engine may be 
significantly different from "open-loop" models of the engine. 

An engine model for the propulsion system is desired which 
describes the engine torque response to fuel inputs from the fuel 
control and to drive train and rotor dynamic loadings (within 
excitable bandwidths). Dynamics in the engine or fuel control 
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which are not excited during this operation are not considered 
important to the modeling task. The sophistication of the engine 
model should be consistent with that of the other elements of the 
coupled rotor/propulsion system model. The coupled model is 
intended to reproduce dynamic phenomena in the zero to once-per- 
rotor-revolution frequency range. 

The value of the engine model will be in its description of 
important aspects of performance. Steady-state values of torque 
and shaft speed should be predicted accurately. Dynamic response 
to small perturbations should accurately reflect engine 
behavior. Match of intermediate variables such as internal 
temperatures and pressures describing engine operation is not 
critical to predicting interaction problems between the rotor and 
propulsion system. 

Models used to describe typical engine behavior vary widely 
in complexity and accuracy [8]. The simplest description of 
engine response is a plot of corrected engine thrust versus power 
lever angle. A typical characteristic is shown in Figure 4.1. 
This plot can be used as a dynamic model if an appropriate time 
lag is associated between the actual throttle position and a 
lagged or virtual position. The engine is observed to accelerate 
and decelerate at rates dependent on the power level. This 
effect can be represented in the model as a variable rate limit. 
The experimental model thus generated is shown in Figure 4.2. No 
attempt is made at a phenomenological explanation of the 
behavior. Parameters are adjusted from observations. Since this 
model does not reflect the internal cause of the torque response, 
the match between the system behavior for various size inputs and 
for different starting and ending values of torque will be poor. 

A more significant drawback of these models is the poor 
closed-loop description of the behavior. This mismatch is 
typically manifested by an incorrect prediction of the 
closed-loop stability of the system. This is due primarily to 
matching the step response in estimating parameters rather than 
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Figure 4.2 Experimental Model of Engine Dynamics 
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matching frequency response over a suitable bandwidth. 

Generally, frequency matching requires more complex model forms. 
In order to accommodate this lack of precision, simple models 
tend to be ’’improved" by ad hoc additions which attempt to 
correct the fundamental non-phenomenological characteristics of 
the formulation. 

Far more complex digital simulations can be generated which 
model component characteristics measured from rig tests and 
aerothermodynamic phenomena occurring in the gas path. Such a 
simulation is shown schematically in Figure 4.3. These 
simulations include basic physical laws relating energy, work, 
mass flow, and acceleration as they dynamically interact in the 
engine. Various "adjustable" parameters such as lumped 
isentropic efficiencies and areas are adjusted to match the 
observed steady-state relationships between input and output 
variables. A detailed representation of the governor is 
included, and the overall response can be tuned to match observed 
dynamics. While the type of simulation recreates steady-state 
performance accurately, the match of transient response is poorer 
due to modeling uncertainties in the complex equations. 

Predicted stability characteristics for closed-loop control are 
quite good, however. 

The drawbacks of detailed analytical simulations are that: 

(1) a large computer capability is necessary; 

(2) detailed internal engine and control information 
is required which is often not available; and 

(3) the added complexity cannot be justified or 
validated from input/output performance observed 
in operation. 

Clearly, a middle ground must exist. 

An experimentally validated analytical description of the 
propulsion system is ideally suited to an analysis of an 
integrated rotorcraft system. It can be developed by considering 
the available measurement data and the type of inputs provided. 
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The development starts with a complete analysis of possible 
dynamic effects. Data are examined and unimportant or unexcited 
phenomena are removed. The remaining model terms are examined 
for bandwidth of response relative to the simulation 
requirement. Frequency effects outside the drive train and rotor 
bandwidths can be neglected. The resulting system represents a 
description of the response including both steady-state and 
dynamic performance whose parameterization is available from 
recorded measurements. This type of model can provide an 
excellent component module in a simulation used to analyze an 
integrated rotorcraft engine, fuel controller, drive train, and 
rotor system. 

4.1.2 Modeling 

The dynamics of the engine in the region near full power are 
nearly linear. Small perturbations in the state of the engine 
are characterized by changes in the speed of the power turbine 
and gas generator turbine. Previous analyses of engine dynamics 
have verified the linearity and the structure of the linear 
system representing this behavior [9-11]. 

Since rotorcraft engines typically have no variable geometry 
in the gas path , the primary dynamic exciter xs the fuel flow 
which is metered by a hydromechanical governor. Fuel is 
chemically converted to heat in the combustor. This provides 
excess energy at the turbine entrance. The combustor lag is 
typically much faster than the dominant system response. The gas 
temperature at the entrance to the turbines is converted into 
work by the expansion through the rotating turbines. This energy 
xs used to drive the compressor on the fixed shaft and the rotor 
on the free shaft. The dominant engine dynamical states can be 
associated with the rotating inertia of these elements. Newton's 
law for the shafts can be written in terms of the speeds and 
torques as follows: 
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where ( )^^ ( 3 ^ refer to the gas generator and power 
turbine turbine shafts. The torques are determined from the mass 
flow and pressure change across the component. The form of the 
expression is given below (for a compressor) 


Tij^m(PRy - 1/Y - 1) 

^COMP " 

where the compression ratio, Pj^, is determined by flow 
equilibrium through the engine and duct. 


(4.3) 


These equations can be used as the basis of a global 
nonlinear model of the engine. Maps of nonlinear parameter 
variations and detailed dynamical representations of 
nonequilibrium flow can model the continuity, energy, and mixing 
phenomena occurring in the gas path. However, these terms can be 
considered negligible if one is interested in system bandwidths 
characteristic of rotor and drive train dynamics. 


For this study, modeling the engine behavior during moderate 
speed excursions by a linear set of dynamical equations provides 
an adequate representation. The form of the system equations can 
be simplified to second-order for most engines by considering the 
physical character of the motion. The linear system of constant 
coefficient equations which is constructed for the engine is 
shown in a Laplace transform block diagram in Figure 4.4 and 
collected in state variable form as follows: 
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Figure 4.4 Linear Engine Model 
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and 


where the five linear parameters, t t Cg„, Cp^y, 

r S Or o 

a, represent the dynamical behavior of the engine. 

The dynamical equations are written for perturbations away 
from an equilibrium condition. A typical choice for this 
equilibrium point is 100 percent power. 

The engine is not the only source of nonlinearity. The 
throttle linkage, electrical, and hydrodynamical fuel controls 
also represent limiting elements in the linearity of the 
response. The components are quite complex. Accurate 
representations of the internal dynamics are difficult to model 
with linear equations. 

The hydromechanical fuel control is the most complex 
nonlinear element of the system. The fuel control meters fuel to 
the engine in response to throttle inputs and engine demand. 
Instability, temperature, speed, and airflow constraints are 
accommodated while producing specified performance at different 
temperatures, altitudes, and speeds. 

Valid dynamic models describing the fuel control hardware 
need to be developed in a parameterized form consistent with Eq. 
(4.4). Yet, because of the control nonlinearity, this is not 
reasonable. The combination of nonlinear control equations with 
Eq. (4.4) will result in an engine/fuel control model suitable 
for analysis of the integrated rotorcraft system for the present 
study . 

4.1.3 Data 


SCT obtained the nonlinear, non-normalized versions of a 
small, 1500-shaft horsepower helicopter turbine engine and fuel 
control system. The simulations used relatively simple engine 
and control models and appeared to provide all the information 
required for the coupled model. 
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The engine simulation used four states: 

(1) power turbine speed, Np; 

(2) gas generator speed, Ng ; 

(3) gas generator turbine inlet pressure, P41; and 

(4) power turbine inlet pressure, P45. 

(Subscripts 41 and 45 on the pressures correspond to stations 
along the engine as were shown in Figure 2.3.) The model used 
nonlinear maps and detailed dynamical representations to 
represent the continuity, energy, and mixing phenomena occurring 
in the turbine gas path. 

The fuel control system consisted of two parts. The first 
part represented the hydromechanical part of the fuel control, 
and the second part represented the electromechanical part of the 
fuel control system. 

The sections that follow describe the parameterization of 
the engine and fuel control models, for the coupled rotorcraft 
model, from this data. 

4.2 ENGINE MODEL 

A four-state linear perturbation engine model was derived in 
the present study to describe the performance of the small 
helicopter engine during moderate speed excursions from an 
operating trim of 100% Np and 97% Ng . This corresponds to 
a 356-ft lb torque output. This model of the engine was 
consistent with system bandwidths characteristic of rotor and 
drive train dynamics. 

The four states used in the model were the same ones used in 
the nonlinear simulation. The nonlinear maps and complex 
dynamical expressions used in the nonlinear model were to be 
replaced by constant coefficients evaluated about the operating 
trim condition. Parameters for the linear model were to be 
extracted from the nonlinear simulation by using the central 
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difference, state perturbation method. The outputs from the 
perturbation model were those required by the engine fuel control 
and by the drive train and fuselage dynamics models. 

4.2.1 Parameterization of Generic Engine Model 

4.2. 1.1 Parameterization Method 


A state perturbation method was used to parameterize the 
linear engine model. The state perturbation approach is 
described in Figure 4.5. 

This method involves freezing the integrators in the 
nonlinear simulation after steady state is reached, then 
perturbing the states and controls one by one to compute the A, 
B, C, and D matrices of the linear state variable representation 
of the system. The approach is shown schematically in Figure 
4.5. 


In Figure 4.5(a), a block diagram of a linear model is 
shown. A linear scalar model is illustrated and discussed, but 
extension to a multivariable nonlinear formulation is 
straightforward . 


The procedure for computing the parameters of the model is 
to integrate the nonlinear model to the steady-state operating 
condition, where x=0. This implies that the system being 
modeled must be a stable, damped system. The integrators are 
then suppressed, as shown in Figure 4.5(b), and a perturbation, 
AX, added to x. The perturbation will cause perturbations Ay 
and AX in y and x. The coefficients in the A and C matrices 
can be computed as 




A similar perturbation, Au, can be input to u to compute B 
and D (see Figure 4.5(c)). The procedure can be extended to 
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A. LINEAR MODEL 


0 AT STEADY STATE 



B. STEADY-STATE. INTEGRATION-SUPPRESSED, PERTURBED x 


. 0 AT STEADY STATE 



C. STEADY-STATE, INTEGRATION-SUPPRESSED, PERTURBED u 


Figure 4.5 Linear Extraction by State Perturbation 
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vector X and u by independently incrementing each element of 
the X and u vectors. 

The magnitude of the perturbations should be consistent with 
amplitudes expected in the real system. Results will generally 
vary with the magnitude of the perturbations used. 

The approach is applicable to stable, damped systems where 
the integrators in the nonlinear simulations of the system are 
accessible and can be suppressed. The procedure may be tedious 
for large order systems. The procedure results in linear, 
dynamic equations in the perturbations 

Ax = A Ax + B Au 

Ay = C Ax + D Au 

4. 2. 1.2 Four-State Engine Model 

A four-state linear perturbation engine model was extracted 
from the nonlinear engine simulation first. Figure 4.6 shows the 
form the linear equations were written in so they could easily be 
connected to fuel control and drive train models. The 
coefficients of the matrices used in the model are given in 
Figure 4.7. 

The operating trim about which this linear engine model was 
developed was 100% Np and 97% Ng at an altitude of 1000 feet 
and an ambient temperature of 34°F. This corresponds to a 
3S6-ft lb torque output. 

Figures 4.8(a) to (f) show the validation of the linear 
engine model with the nonlinear engine model results. As can be 
seen the fit is very good up until the engine reaches 110% shaft 
speed at time 15 seconds. The open-loop eigensystem of this 
linear engine model is shown in Table 4.1. All four eigenvalues 
of the model are non-oscillatory , exponentially decaying modes. 
The time constants of the modes range from 12.4 times the rotor 
period to 0.02 times the rotor period. The fastest mode ( '^ = 
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WHERE 


AXg = A^£ AX£ + 

Ay£ = C££ Ax£ + DgpAyp 

^E * {%r ^g* '’as* 

= ENGINE STATES 

^E * {^E * ''’45’ ^S 3 * ’^ASr} ^ 

= ENGINE OUTPUTS 

yj » {op^, Qp2* ^tR» ^MR^ 

= DRIVE TRAIN OUTPUTS 

= {Wf} 

= FUEL CONTROL OUTPUTS 

Aee» B£j, Bgp, C££, D£p = MATRICES OF COEFFICIENTS 
A( ) = PERTURBATION VALUE 


Figure 4.6 Four-State Engine Model 
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Figure 4.7 Coefficients for Four-State Linear Perturbation 

Engine Model 
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Figure 4.8a Validation of Four-State Linear Engine Model 
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Figure 4.8b Validation of Four-State Linear Engine Model 
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Figure 4.8f Validation of Four-State Linear Engine Model 






Table 4.1 

Eigensystem o£ Four-State Engine Model 


NO. 

EIGENVALUE 

TIME CONSTANT 
(SEC) 

MODE DESCRIPTION 

1 

-.28 + Oi 

3.6 

Power Turbine Speed (Np) 

2 

-1.4 + Oi 

.71 

First Turbine Stage Inlet 
Pressure (P 41 ) 

3 

1 

• 

00 

+ 

o 

.13 

Gas Generator Speed (Ng) 

4 

-150. + Oi 

.007 

Power Turbine Inlet Pressure 

(P 45 ) 


NOTE: ^ = .028 SEC/CYCLE 

“mr 
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.007 sec) IS faster than the dynamics required of the generic 
rotorcraft model and was determined to be unneccessary . 

4. 2. 1.3 Model Reduction 

Model reduction was applied to the four-state linear engine 
model to simplify it. Two states were proposed (see Section 
4.1.2, Eq. (4.4)) for the generic model structure, but four 
states were used originally. It was reasonable that the number 
of states in the four-state model be reduced since the model 
contained higher frequency characteristics than necessary, and 
would be inconsistent with modeling assumptions of the other 
subsystems. By setting the derivatives of faster states to zero, 
model order reduction was used to remove the faster mode. This 
approach is described in Figure 4.9. 

Figure 4.9 illustrates model order reduction for a two-state 
model. The first equation in the figure is the state equation 
for the system. Suppose it is found analytically (or in an 

eigenanalysis) that X 2 is much faster than X]^ (i.e. 

• • 

^2 Then X£ is assumed to change instantaneously 

with respect to x^, and X£ is set to zero. The equation 
for X 2 can then be solved for X 2 in terms of x^ and 
u. This result is substituted into the equation for x^ to 
remove X 2 from the equation (as shown in the last equation of 
Figure 4.9). The approach can easily be extended to a many 
degree-of -f reedom system and also applied to remove unwanted 
state(s) from output equations. 

Table 4.2 shows the eigensystem of the three-state model 
that resulted after model order reduction was applied. Note that 
the eigensystem for the three states left in the model was 
practically the same as that of the three slowest modes before 
the fast mode was removed (see Table 4.1). Plots of the 
responses of the three-state model were indistinguishable from 
those of the four-state model for the inputs and plotting scales 
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WHERE: MUCH FASTER THAN 

SET x^ = 0 

^2 = 0 = [321 ^2 M " 

OR 

x^ = -322 32^ - 322 b2 ^ 

THUS 

^1 *-®lT®12 ^22 ®21^ l^l' ^ ^*^1“ *^2^22 ^2^{'^} 


Figure 4.9 Model Order Reduction 
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Table 4.2 

Eigensystem of Three-State Engine Model 


NO. 

EIGENVALUE 

TIME CONSTANT 
(SEC) 

MODE DESCRIPTION 

1 

-.28 + Oi 

3.6 

Power Turbine Speed (Np) 

2 

-1.3 + Oi 

.77 

First Turbine Stage Inlet 
Pressure (P 41 ) 

3 

-8.0 + Oi 

.13 

Gas Generator Speed (Ng) 


NOTE: ^ = .028 SEC/CYCLE 
“MR 
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used in Figures 4.8(a) through 4.8(f), so no new plots for the 
three-state model are shown. 

Figure 4.10 shows the form of the generic model structure 
for the three-state engine model. Coefficients for the model are 
given in Figure 4.11. 

It is reasonable that only one pressure state be used in the 
three-state model. Theoretical fluid mechanics indicate that a 
quasi-steady relationship exists between pressures. Any 
disturbance of the relationship among pressures is quickly 
corrected, so only one pressure state is needed. 

4. 2. 1.4 Engine Subroutine 

The three-state engine model become one subsystem component 
of the rotorcraf t/propulsion simulation. Trim values were added 
back onto the perturbational engine states and output variables 
in the model so total as well as perturbational outputs would be 
available. Figure 4.12 shows a flowchart of this engine 
simulation. 

4.3 FUEL CONTROL MODEL 

The fuel control simulation with the nonlinear engine 
simulation was converted directly into a subroutine for use in 
the generic rotorcraf t/propulsion model. It was not linearized 
or simplified. A flowchart of the fuel control model subroutine 
is shown in Figure 4.13. This fuel control simulation could be 
used to control either one of the linear or the nonlinear engine 
simulations . 

4.4 DRIVE TRAIN MODEL 

Figure 4.14 shows the proposed generic structure for the 
drive train model. The equations for this model are described in 
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WHERE 


A^e = Agg Ax^ + Bg/yy + B^p Ayp 
AyE = Cpp Axp + DppAyp 


^E “{*’41’ '^g’ 

= ENGINE STATES 

^E " {^E ' ''’45’ '’SS’ '^45R I"*" 

= ENGINE OUTPUTS 

^T {^Pr ^P2’ ^TR* ^Mr} ^ 

= DRIVE TRAIN OUTPUTS 

yp = { wp} 

= PUEL CONTROL OUTPUTS 

%* Bet. ®ep‘ ^ee* °ep " matrices op coeppicients 

A{ ) = PERTURBATION VALUE 


Figure 4.10 Three-State Engine Model 
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Aee = 


Bet = 


Bef = 


Cee - 


Def = 


Figure 


-3.95 .0425 0. 

250. -5.36 0. 

L44.7 -.031 -.283 J 


0 0 0 

0 0 0 

L-15.92 0 0 


0 

0 

0 


'1260. 

25500. 

L 12100. J 


2.8 -.00196 -.0178 I 

.187 -.0432 0. 

.926 .00156 0. 

.000448 -2.35E-6 0. 


r 762. T 
5590. 

0 . 

L .385 


.11 Coefficients for Three-State Linear Perturb- 
ation Engine Model 


58 




Figure 4.12 Flow Chart of Generic Engine Subroutine 
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Figure 4.13 Flow Chart of Generic Fuel Control Subroutine 
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WHERE aX^ 


“mr| 

> => DRIVE TRAIN STATES 


AY 


T 



DRIVE TRAIN OUTPUTS 


AXg^ . aX^2 ’ engine STATES FOR TWO ENGINES 
AY^ - ROTOR OUTPUTS 
AY|_« LAG DAMPER OUTPUTS 

*TET Ve2’ *TT’ ®TR’ ®TL’ ''TT’ ^TET ^TE2 ' COEFFICIENTS 


Figure 4.14 


Form o£ Generic Drive Train Model 
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Figure 4.15. The derivation of the equations is straight forward 
and similar to that used to derive the equations in [12]. 

The drive train model is shown schematically in Figure 
4.16. Note that the transmission inertia has been lumped into 
the rotor hub inertia. This simplification was made because 
numerical data were not available for modeling the transmission 
separately. The power turbine inertias shown in the figure are 
not really in the transmission equations, Figure 4.15, because 
they were already included in the engine model equations. 

The parameters used in the model are representative of a 
small, twin turbine helicopter. They do not correspond to any 
one helicopter. Although this drive train is not the one 
connected to the engine, fuel control, and rotor being used in 
consideration, it is thought to be sufficiently representative 
for the purposes of demonstrating the simulation model. Figure 
4.17 lists the parameters used in the model. 

Damping parameters were not available. Small values of 
damping were used to assure that the system would be stable. 

4.5 MODEL CHECKOUT 

The propulsion system model was checked out by executing it 
independent of the rest of the rotorcraft simulation. The 
propulsion system that was constructed includes two engine 
subroutines, two fuel controllers, and a drive train subroutine. 

A rotor load torque time history was specified to excite the 
propulsion system dynamics. The results of a checkout run are 
plotted in Figure 4.18. The plots in this figure show the 
responses of the propulsion system with both full nonlinear and 
three-state linear engine models. 

It might be useful in the future to perform an eigenanalysis 
of the engine/drive train/engine part of the simulation to 
determine its frequency content. It might be shown possible to 
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STATE 



Figure 4.15 Equations For Drive Train Model 




Figure 4.16 


Drive Train Model 


dg =0. 

dH = 0. 

dnE = 1070. ft lb sec (referenced to hub) 
Ih = 880. ft lb sec2 (referenced to hub) 
kwE = 64930. ft Ib/rad (referenced to hub) 
rm = TH2 = .010572 


Figure 4.17 Drive Train Model Parameters 
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Figure 4.18c Linear and Nonlinear Propulsion Systems 
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reduce the order of this system by a few states, and still cover 
the desired frequency range. 

Additional system dynamics, not demonstrated in the plots of 
Figure 4.18, are the engine hunting modes. Future propulsion 
system runs should be performed that perturb one engine 
differently from the other, and thus excite engine hunting. The 
damping of the engine hunting mode could be observed after the 
excitation was removed. 

4 . 6 SUMMARY 

A rotorcraft propulsion system model of a particular engine, 
fuel control, and drive train combination was constructed. A 
methodology for parameterizing linear versions of the model was 
discussed, and parameters for a linear engine model developed. 
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V. ROTOR SYSTEM MODELING 


This chapter describes a methodology for constructing a 
rotor/fuselage model. The form of the generic model structure is 
also described. The model is based on the structure defined in 
Chapter III. Model parameters of a particular rotorcraft, used 
to demonstrate the model, are obtained from a nonlinear 
simulation of a research helicopter. Regression techniques are 
used to obtain these parameters. 

5.1 ROTOR/FUSELAGE GENERIC MODEL FORM 

The rotor/fuselage generic model equations were constructed 
in the input-output form shown in Figure 3.3. This form allowed 
the model to be plugged into the modular framework for the 
coupled system that was described in Figure 3.2. This 
input-output approach also allowed the rotor and fuselage models 
to be developed separately from the other modules (e.g. engine, 
drive train) of the rotorcraft system. 

The states used in the rotor model were multi-blade 
coordinates (e.g. Bq , Cq, ?^(-,...). The states 

in the fuselage model were the six rigid-body rates (u, v, w, p, 
q, r) and the three rigid-body Euler angles ( <|), 9, 'P) . 

The rotor and fuselage models were interconnected through 
the rotor hub forces and the fuselage states. The forces were 
outputs of the rotor model. 

The forces (and moments) acting on a helicopter fuselage 
include aerodynamic forces, main rotor forces and tail rotor 
forces. These forces were kept separate in the model to provide 
a better framework for identifying the structure of these forces 
and to provide flexibility in the model. 

The dynamics of the fuselage and rotor model and the 
fuselage aerodynamic forces were described using perturbational 
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quantities, as was done for the propulsion system. Trim values 
were added back onto the perturbational quantities in the 
subroutines, as was also done for the propulsion system. The 
trim values were added back on after integration so total 
response will be available. 

5.2 GENHEL RSRA SIMULATION 

The GENHEL RSRA simulation was used in parameterizing the 
perturbational rotor and fuselage models. It is a fully 
nonlinear simulation of the main rotor, tail rotor, fuselage and 
control system of the Rotor Systems Research Aircraft (RSRA). It 
does not include rotor propulsion system dynamics (i.e. engine, 
fuel control, drive train dynamics). The equations in the 
simulation originated from Sikorsky Aircraft's General Helicopter 
Simulation (GENHEL) which is documented in [13,14]. The main 
rotor mathematical model in the simulation is a rotating blade 
element model of an articulated rotor system. Rotor blade and 
airframe flexibility are not taken into account. Flapping and 
lagging hinges are assumed coincident. A Bailey model [15] is 
used for the tail rotor. 

The control system represented in the simulation includes a 
full complement of helicopter controls. It includes the pilot's 
control system and the stability augmentation system (SAS). 

The GENHEL simulation has a modular format, as shown in 
Figure 5.1. It is well documented and is easily accessible on a 
CDC-7600 batch computer. 

5.3 PARAMETERIZATION METHODOLOGY 

The parameters (or coefficients) for the rotor and fuselage 
equations were obtained from time history data generated by the 
nonlinear GENHEL simulation. A stepwise regression algorithm, in 
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use at SCT, called OSR was used to compute the coefficient 
values . 

The stepwise regression method, OSR (Optimal Subset 
Regression), requires as input the time histories of the states, 
state derivatives, and control inputs. These inputs are the data 
that can be generated by the GENHEL simulation. 

Conceptually, the use of OSR involves setting up an equation 
for each acceleration term, mx, as 

mXi = EF^ 

= Cl Xi + C^X^ * C3U + . . . 

where the terras x^, xj^, u, etc., are candidate independent 

variables and the coefficients C2» Cj...) are 

•• • 

unknown. The variables x^, xj, x^, U]^, etc., are input 
to OSR as time histories. OSR then determines the subsets of the 
candidate independent variables which best represent the data for 
the dependent variable, mx^. OSR also computes estimates of 
the unknown coefficients of the independent variables in the 
subset. The subset independent variables and variable 
coefficients are chosen to minimize the mean square error between 
the mX]^, time history and the subset time history. Thus, the 
"fit" of the linear equations is a mean square minimization. 

The stepwise regression procedure implemented in OSR is not 
applied to data from the steady state performance of the system, 
but is applied to data from specific trajectories or maneuvers. 
The trajectories are defined by sets of control inputs chosen to 
excite the system performance that is to be modeled. 

A single maneuver is often insufficient to properly excite 
all the states so a sequence of maneuvers is performed with the 
rotorcraft returned to its trim condition between maneuvers. 

Subset regression has an advantage over perturbation methods 
in that it automatically averages over a time interval or 
period. This averaging is important when identifying constant 
coefficients for a system that has significant periodic 
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coefficient effects. An additional advantage of the subset 
regression method is that it can easily handle a large number of 
states . 

The stepwise regression method as formulated above applies 
to a linear or non-linear model of the motions of a dynamic 
system from test or simulation data. The approach can also be 
used to determine a perturbation model (i.e., where the states 
are perturbation motions about a nominal or trim state rather 
that total motions). 

Difficulties are often encountered in the application of 
this approach when there is a high correlation between any of the 
states and/or controls. Correlation causes a reduction in the 
linear independence of the regression variables and degrades the 
identification process. 

5.4. ROTOR MULTIBLADE EQUATIONS 

The perturbational rotor model was defined in the 
multi-blade coordinate system. Recall that the nonlinear GENHEL 
simulation defined the rotor blade dynamics in the rotating 
coordinate system. Multiblade coordinates are preferred for the 
perturbational model because they are constant in steady flight 
(except for higher harmonics) and because they provide accurate 
information on the magnitude and direction of the thrust of the 
rotor. They are useful in modeling and gaining a better physical 
understanding of the coupling between the rotor and the 
fuselage. They are also useful in feedback control synthesis. 
Rotating blade coordinates on the other hand allow a better 
physical understanding of the forces on, and dynamics of, a 
single blade, but do not provide for easy understanding of the 
effects on the fuselage. 

The transformation from the rotating blade coordinate system 
to the multi-blade coordinate system is described in Appendix A 
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for a five-bladed rotor system. This appendix also discusses the 
form of the resulting equations which were parameterized. 

5.5 ROTOR/FUSELAGE EQUATIONS 

The final forms of the fuselage aerodynamic and rotor 
equations are shown in Figures 5.2 and 5.3. More detailed 
descriptions of the matrices in the rotor equations are given in 
Appendix A where the multi-blade transformation is discussed. 
Parameterization for the models is discussed in the following 
section. 

The unknown or uncertain portion of the fuselage dynamics is 
the aerodynamics. Rigid body dynamics are well understood. 
Standard simulations are available. A regression to identify a 
fuselage model need only identify the particular form of the 
aerodynamic equations. 

In the current study the fuselage aerodynamic forces are 
separated from the hub forces and identified as a separate module. 

5.6 PARAMETERIZATION OF ROTOR AND FUSELAGE MODELS 

Parameters for the linear rotor and fuselage aerodynamics 
models were obtained by applying OSR to the perturbation data 
generated by executing the GENHEL simulation. This section 
describes this work. 

The trim condition about which the GENHEL-RSRA simulation 
was operated was as follows: 

• forward speed = 120 knots 

• altitude = 1000 feet 

• gross vehicle weight = 19,600 pounds 

• ambient temperature = 59°F 

• ambient pressure = 14.7 psi 
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AX(^ = Aj^(^ AXpj + AX-J. + Aj^p^j AXpjj + AU + Bj^py Aypy 


^ Srt ^ Srl ^ 


,(J - J'^) 2J J 


AX, 


‘J'R = "^RR “r * '=RT ‘*T * ‘^RFU “'’'FU * "rT * °RC 
* ''rRI ''^IR * ’’rFU *^FU 


where 

\r’ ^RT’ Vu’ ®RC’ ®RFU* ®RT’ ®RL’ •'RR* ^RT’ ^RFU’ 

^RC’ ^RRl* ^RFU “ matrices of coefficients 

J = axis transform matrix (see Eq. (A.4) in Appendix A) 

u = rotor controls 

= i&c. Axs, Bis. »TR ^ ^ 

xpu = fuselage states 

» |u, V, w, p, q, r, d, s, 

XR = rotor states 

~ I ®o ♦ ***’ ®2s> ^^o> •••» ^2s» 8o» •••» 

®2s» Co> •••1 ^s 

xj = drive train states 

yFU * fuselage outputs 

= |u, V, w, p*, q, r, /v/ } T 

yR = rotor force and torque outputs 

= I H, J, T, Lh, Mh, Qh-QlD f ^ 

yj = drive train outputs 

yiR ~ }®o» •••» B2s>Co» •••> C2sf"^ 


Figure 5.2 Form of Rotor Perturbational Equations 
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’ '•AFU *’'fu * Sr ‘''‘r * Sc * ‘’at *■>'1 * Sfu ‘'^fu 
■" “aRI ‘^IR 

where: 

^AFU* ^AR* °AC’ “^AT’ °AFU> °AR1 = ^oeff^c^ent matrices 


u 

3i 

rotor controls 

XFU 

= 

fuselage states 

|u, V, w, p, q, r, 4>, e, [T 

XR 

= 

rotor states 

yA 

S3 

fuselage aerodynamic force outputs 
|xF. yF. ZF, Lp, Mp, Np fT 

yFu 

a 

a 

fuselage outputs 

|u, V, w, p, q, r, /v/|T 

yj 

a 

drive train outputs 

yiR 

a 

rotor accelerations 


Figure 5.3 Form of Fuselage Aerodynamicss Perturbation 
Equations 
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The data generated by the GENHEL simulation was stored on a 
computer disk to be used in the regression analysis. These data 
included the states, state derivatives and control inputs. 

5.6.1 Control Inputs 

The control inputs were defined so they would excite all of 
the rotor and fuselage dynamics in the frequency range 0 to 
2 where is the main rotor speed. Five different 

control inputs were used. These included 

- main rotor collective 

- main rotor lateral cyclic 

- main rotor longitudinal cyclic 

®TR ■ rotor collective 

Qg - engine torque on rotor hub 

The engine torque, Qg, on the rotor hub, was included 
because the GENHEL simulation did not have a rotor propulsion 
system. The engine torque was used to vary the rotor speed and 
to thus excite the rotor dynamics that needed to couple with 
propulsion system dynamics. 

A doublet was used as the form for the control input 
perturbations. Doublets at four different frequencies were used 
in order to excite all response frequencies of the rotor and 
fuselage between zero and The shape and duration of 

these doublets are shown in Figure 5.4. 

5.6.2 Attitude Control 


The GENHEL simulation contained a Stability Augmentation 
System (SAS) to stabilize the fuselage dynamics. The SAS would 
not "fly" the simulation, though. Some type of pilot-in-the-loop 
was required to keep the helicopter from flying off erratically. 

A simple attitude controller was .coupled with the SAS in order to 


79 



(a) TRIM CONDITION 

• ALTITUDE = 1000 FEET 

• FORWARD SPEED = 120 KNOTS 
t TEMPERATURE = 59°F 
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(b) INPUT FORM 
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Figure 5.4 Generation of GENHEL Data for OSR 



keep the helicopter near the nominal trim flight condition that 
was being modeled. 

The attitude controller consisted of a simple proportional 
controller applied to the roll and pitch axes of the fuselage. A 
schematic of the controller for the roll axis is shown in Figure 
5.5. The gains for the roll and pitch controllers were chosen to 
be small so as not to wash out the transient effects of the 
control perturbations. 

5.6.3 Data Generation 


Table 5.1 summarizes the data runs used to generate data for 
the regression analysis. Twenty data generation runs were made 
(five different controls times four doublet periods). The run 
times and sampling rates are also given in Table 5.1. The data 
from these twenty runs was catenated end-to-end to create one 
long time history of 3950 time points. 


5.6.4 Dependent and Independent Variables 


Table 5.2 summarizes the dependent and independent variables 
that were used in the regression analysis. The independent 
variables are the ones that are collected into equations to 
describe each dependent variables time history. Table 5.3 
describes the nomenclature used in Table 5.2. 

Some additional constraints, pertaining to which dependent 
variables were allowed into the equations, were applied that are 
not shown in this table. These constraints forced the diagonal 
matrix form, found in Eq. (A. 10) of Appendix A, into the 
regression equations. 
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Figure 5.5 Attitude Control For Roll Axis 
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Table 5.1 

Description o£ Data Runs Used to Generate Data for OSR 


a. CONTROL DISTURBANCES 


NO. 

CONTROL 

NOMINAL VALUE 

PERTURBATION 

1. 

©C- COLLECTIVE 

11* 

1.5* 

2. 

Ais- LAT. CYC. 

-1.5* 

0.3* 

3. 

Bis- long. CYC. 

9.3* 

0.9* 

4. 

©TR- TAIL ROTOR COLL. 

5.9* 

-2.2’ 

5. 

Qe- eng. torq. on hub 

31630 FT. LB. 

3160 FT. LB. 


b. DISTURBANCE FORMS 


NO. 

^1 

m 

SAMPLE INTERVAL 
T (SEC) 

1. 

.142 ( 2/n) 

.5 

.02 

2. 

.284 { 1/n) 

1.0 

.02 

3. 

1.14 (.25/n) 

6.0 

,05 

4. 

5.68 (.02/n) 

20. 

.10 
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Table 5.2 

Regression Model Setup 


DEPENDENT VARIABLE 

INDEPENDENT VARIABLES 

AXp 

- GROUP 1 (au, av, aw, ap, Aq, at) 

aYf 


aZf 

- GROUP 2 (aSq,..., a62s» a5o.->*. 4C2$t 

aLf 

• * * * \ 
ABq,..., A82s» ^^o»***» ^^2s’ 

aMf 

GROUP 3 (^6q9*«*9 a02s* ^^0***** 

aNf 

- GROUP 4 (au, av, aw, ap, aq, at, a/v/) 


- GROUP 5 (a^tr, AfJjR) 

t> 

: 

o 

- GROUP 1 

*ABlc 

- GROUP 2 

*A6ls 

- GROUP 4 

*AB2c 

— GROUP 6 (a6^, aA^5, aB^s* aS^mr) 

*A 82 s 


*ACo 


*Ak’lc 


*A?ls 


*4C2c 


*AC2s 


AH 

- GROUP 1 

aJ 

- GROUP 2 

aT 

- GROUP 4 

aLh 

- GROUP 6 

aMh 

- GROUP 7 (a^mr) 

aQh - 4Qld 



84 









Table 5.3 


Nomenclature For Table 5.2 


A ( ) - PERTURBATION VALUE 

*( ) - DOES NOT INCLUDE AXIS TRANSFORMATION "J" TERMS 
{SEE APPENDIX A) OR LAG DAMPER TORQUES 

Xp,..., Np . AERODYNAMIC BODY-FIXED FORCES AND MOMENTS 

ON FUSELAGE 

®0»---®2S - flapping MULTIBLADE COORDINATES FOR 5 BLADES 

" RO'J’OR lagging multiblade coordinates for 5 BLADES 
H Qh - HUB FORCES ON ROTOR IN SHAFT COORDINATES 

Old - lag damper torque for 5 blades 

u, V, w, p, q, r -FUSELAGE RATES IN BODY-FIXED COORDINATES 

®C. Aic, Bis - MAIN ROTOR COLLECTIVE, LATERAL CYCLIC AND 

LONGITUDINAL CYCLIC CONTROL INPUTS 

®tr - tail rotor collective input 

^MR. ^TR - main rotor AND TAIL ROTOR SPEEDS 
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5.6.5 Results 


5.6. 5.1 Coefficients 

The coefficients for the fuselage and rotor models that were 
obtained from the data are listed in Volume II. 

Identification of the fuselage aerodynamic models (equations 
for Xp.Yp,.., Np) and for the rotor hub forces and 
moments (equations for H, J, T, L^, Q^) were 

straightforward and proceeded without difficulty. OSR runs were 
made using the dependent and independent variables listed in 
Table 5.2. Very good fits to the data were obtained. 

Significant difficulty was encountered in attempting to 
identify coefficients for the rotor dynamics model. Initial OSR 
runs were made with the "axis transformation" terms (the terms 
involving "J" in Eq. (A. 8) of Appendix A) removed from the 
regression. "Removed from the regression" means the terms were 
known beforehand and were subtracted off of the acceleration term 
before attempting to fit an equation to the acceleration term. 

All of the blade states, fuselage states, controls (except for 
the tail rotor), and fuselage accelerations were allowed in the 
analysis during the initial runs. No attempt was made to force 
the diagonal matrix coefficient form of Eq. (A. 10) of Appendix A 
into the analysis. The analysis results of Appendix A were not 
available until the completion of this study. 

The model identified in the initial regression analysis 
(just described), had very high mean squared fit values but was 
unstable when integrated. Since the nonlinear model equations 
were stable, the identified model was definitely wrong. Although 
the mean squared (and also the F-test) fit measures said the fit 
was good, it was not. It is known, though, that when there is no 
uncorrelated noise in the data, as was true for this data, the 
F-test is of questionable value. Also, the mean squared fit 
value increases with each variable allowed into the regression. 
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A good fit, although not necessarily physically reasonable, can 
always be obtained by allowing enough variables into the 
regression . 

Several attempts were made to obtain better equations. In 
the first attempts, the number of blade states was reduced. The 
2c and 2s blade lagging states were dropped, then the 2c 
and 2s blade flapping states were dropped. It was thought that 
perhaps there was insufficient information on these states in the 
data or that identification of these states was corrupted by high 
frequency (greater than responses in the data. No 

significant improvement in the resulting model resulted from 
these changes. 

Other changes in the regression that also yielded no improvement 
were ; 

(1) filtering the data to remove responses with 
frequencies higher than approximately 

(2) using larger control input perturbations to try to 
excite more inplane (lag) response, and 

(3) not allowing the fuselage accelerations into the 
analysis. 

The results shown in Volume II were finally computed by 
forcing the diagonal form of the coefficient matrices, shown in 
Eqs. (A. 8) and (A. 10) of Appendix A, into the regression 
analysis while also analytically subtracting off the fuselage 
acceleration terms before performing the regression analysis. 
These acceleration terms were found analytically by linearizing 
and partially transforming the GENHEL equations into multi-blade 
coordinates. The rotor equations thus constructed using the 
regression analysis were stable, although when integrated, did 
not fit the nonlinear response particularly well. A sample of 
the results is shown and discussed in the next section. 

The reason for the difficulty in obtaining the parameters 
for the rotor equations is now understood. It was uncovered 
during the analysis of the multi-blade transformation discussed 
in Appendix A. 
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The difficulty was caused by periodic effects in the rotor 
dynamics. Generally, periodic effects in the rotor dynamics are 
considered to be small and can be neglected up to rotor advance 
ratios of = 0.3 [18]. The flight condition analysed in this 
contract was a helicopter forward speed of 120 knots, which 
corresponded to an advance ratio of 0.3. Thus, periodic effects 
should have been expected to be significant. 

The problem was that the periodic effects were not handled 
correctly in the initial analysis. As shown in Eq. (A. 11) of 
Appendix A, trim for helicopters at forward speeds where periodic 
effects are present is not time-independent. The trim state is 
given by the sura of a steady, non-accelerating term and a 
periodic, time-varying term. When computing perturbations from a 
trim state, the periodic portion of the trim must be removed as 
well as the time- independent portion. This periodic portion of 
the trim was not removed when the perturbat lonal data for the 
regression analysis was saved. 

The procedure used in this study to subtract off the trim 
values was as follows. First the GENHEL simulation was trimmed 
out. Then second, a "snapshot" of the states, responses and 
controls at this trim condition was taken and stored. Third, 
this "snapshot" of the data was subtracted off of the data as it 
was stored during the helicopter maneuvers to create the 
perturbation data. 

This approach was in error because the "snapshot" of the 
data was distorted by the time-dependent trim effects. This is 
shown in Figure 5.6. The correct procedure would have been to 
subtract off the steady part of the trim value, labelled in 
Figure 5.6, and the periodic trim signal. The "snapshot" could 
occur during any phase of the periodic motion so it could error 
from the steady trim by up to * the amplitude of the periodic 
motion if the trim values are not removed correctly. The 
magnitudes of these amplitudes are observable in the plots in the 
next section. If it is only possible to remove the steady 
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Figure 5.6 Example of Trim in a Periodic System 
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portion of the trim, OSR may successfully average out the rest, 
assuming the steady portion is removed correctly. 

The error in removing the steady trim caused the 
perturbation data that was stored to be biased. Each channel (or 
variable) in the data had a different bias value. 

The correct way to compute the perturbation values would 
have been to save the trim response for a period of time. Then, 
average the data to find the steady trim value and fit a sine to 
the periodic portion of the trim signal. Finally, subtract off 
the steady trim value and the sine from the maneuver data to 
leave the perturbation values. 

There were insufficient resources left in the contract at 
the end to redo the rotor model parameterization and correctly 
consider the periodic effects. This was because the regression 
data had been originally taken incorrectly. 

5.6. 5. 2 Simulation Results 


The linear rotor and fuselage models that were extracted 
from the GENHEL simulation were appended back to the GENHEL 
simulation for checkout. A nonlinear lag damper model in 
multi“blade coordinates was added to complete the rotor model. 

The parameter values in the linear perturbation models are given 
in Volume 1 1 . 

The plots in Figure 5.7 show the comparison of the linear 
perturbation rotor and fuselage results with those of the GENHEL 
simulation. Both sets of results were computed assuming a 
constant rotor speed Ci.e, no propulsion dynamics). The trim 
condition is 120 knots forward speed, 1000 ft altitude, which is 
the same trim for which the perturbat ional data were developed. 
The plots in Figure 5.7 show the response to a pilot longitudinal 
cyclic input. The input is ramped up from 78% pilot stick 
position to 88% over one second. 
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Figure 5.7j Rotor/Fuselage With Long Cyc Input 








The effect of the poor rotor model parameters is apparent. 
Figure 5.7(a) shows the pilots longitudinal cyclic input. Figure 
5.7(b) shows the total longitudinal cyclic input at the rotor 
hub, including pilot input, SAS input and attitude controller 
input. Figures 5.7(c) to 5.7(f) show__some of the blade 
accelerations. The ones shown are *^2 c *^o • 

The symbol implys that the terms are not total 

accelerations, and the "axis transform" or "J" terms in Eq. 

(A. 8) of Appendix A have not been added on. This was the 
approach taken in Ref. 20. 

The periodic effects are quite apparent for the nonlinear 
rotor simulation in Figures 5.7(c) to 5.7(f). These effects have 
a major effect at the = 17.6 Hz frequency. Note that 

time scales on some of the plots have been expanded to allow this 
frequency to be eyeballed. 

The effect of the periodic motion on the choice of the trim 
condition is especially apparent in Figure 5.7(e). Here, trim 
data was taken near the peak of the periodic motion of this 
response . 

The remaining figures, Figure 5.7(g) to 5.7(j) show some of 
the forces and fuselage motions. The periodic effects are mostly 
filtered out before they reach these responses, although they are 
still quite apparent in the rotor hub forces. These four plots 
also show that the perturbat ional model is drifting off in error 
due to the poor rotor fit as time increases. 

5.7 SUMMARY AND RECOMMENDATIONS 

A generic methodology was defined for developing a 
simplified rotor/fuselage model from detailed nonlinear 
simulations. The approach is generic in that it can be applied 
to any helicopter or rotorcraft. The modular structure used in 
the approach also allows parts (i.e. rotor, fuselage or lag 
damper) of the total model to be obtained from different data 
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sources, then combined into a single simulation. The simulation 
itself, used here as a test case, was not generic as it applied 
to one particular twin engine, single main rotor rotorcraft 
simulation. 

There were problems in the parameterization of the rotor 
model. The rotor model did not fit the source data very well. 

The problem was attributed to initial lack of knowledge about the 
effect of periodic rotor dynamics. The periodic effects were 
subsequently studied and found to cause a component of the trim 
solution to be periodic. This periodic trim was not correctly 
accounted for when the perturbation data was taken. 

An analysis of the transformation from rotating blade 
coordinates to multi-blade coordinates was performed. This 
analysis showed the form of the rotor equations in multi-blade 
coordinates and also the existence and equation form of the 
periodic terms. The analysis is documented in Appendix A. 
Periodic effects cannot be modeled by constant coefficient linear 
equations. The new information on the periodic effects can be 
used to obtain better multi-blade coordinate rotor models. The 
periodic effects can be either averaged out, now that their 
structure is known, or included in the model. 
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VI. SYSTEM MODEL FORMULATION 


The rotorcraft model developed in this contract was appended 
to the nonlinear GENHEL simulation to form a single simulation. 
This was done for two reasons. First, the perturbation model 
developed in this contract needed a flight control system to 
"fly" it. The GENHEL simulation had a flight control system 
which could be used for this purpose. Second, combining the 
simulations together facilitated their comparison. Checkout of a 
parameterization of the generic model structure could be 
accomplished by running the nonlinear and the perturbation models 
with the same controller and then comparing the output directly. 

Formulation of the combined perturbation and GENHEL 
simulations into one simulation was divided into two steps. Step 
one consisted of appending nonlinear engine, fuel control and 
drive train subroutines onto the GENHEL simulation. This step 
was taken to provide a full nonlinear helicopter model capable of 
simulating the dynamic interactions between the rotor and 
propulsion systems. The second step was to append the 
perturbation model equations to the GENHEL simulation. 

The sections that follow in this chapter describe (1) the 
addition of the nonlinear propulsion system dynamics simulation 
to the GENHEL simulation, (2) the formulation of the combined 
perturbation and GENHEL simulations, and (3) demonstration and 
discussion of the resulting model. User's notes and common block 
descriptions for the augmented GENHEL simulation are given in 
Volume II. 

6.1 ADDITION OF NONLINEAR PROPULSION DYNAMICS TO GENHEL 

SIMULATION 

The nonlinear engine, fuel control and drive train 
subroutines discussed in Chapter IV were appended to the GENHEL 
simulation on the CDC 7600 at NASA-Ames. This addition created a 
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helicopter simulation that could model the interactions between 
the propulsion system and the rotor system. 

The nonlinear propulsion subroutines (engines, fuel 
controls, drive train) were connected to the ROTOR subroutine in 
the GENHEL simulation. A simplified flow chart of the GENHEL 
simulation was shown in Figure S.l. The propulsion subroutines 
contained their own fourth order Runge Kutta integration 
package. Options were written into the GENHEL simulation so the 
simulation could be run assuming the following three propulsion 
system models: 

(1) constant rotor speed, 

(2) constant engine torque on rotor hub, and 

(3) nonlinear engine, fuel control and transmission 
dynamics . 

Figures 6.1(a) through (i) demonstrate the response of the 
resulting augmented simulation. The change in rotor collective 
due to a ramp input of 20% to the pilots collective control is 
shown in Figure 6.1(a). Figure 6.1(b) shows the rotor speed 
droop. The two curves in this figure correspond, in one case, to 
a constant speed rotor model, and in the other case, to the 
nonlinear engine/fuel controller model with the collective input 
fed forward to the fuel control. The plots that follow these two 
show the fuel flow rate, rotor torque, rotor moments, and 
fuselage motions for the two configurations. The addition of 
propulsion system dynamics had little effect on the response of 
the system except for the effects on yaw rates, as shown in 
Figure 6.1(g). 

A paper [21] written in conjunction with this report 
provides additional results and a discussion of the influence of 
the propulsion system dynamics on this model. 
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Figure 6.1b Nonlinear Rotor/Propulsion System With Collective 

Input, Rotor Speed 




Figure 6.1c Nonlinear Rotor/Propulsion System With Collective 

Input, Fuel Flow 




MAIN ROTOR TORQUE (FT LB) 



Figure 6. Id Nonlinear Rotor/Propulsion System With Collective 

Input, Rotor Torque 
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Figure 6.1e Nonlinear Rotor/Propulsion System With Collective 

Input, Body Axis Pitch Moment 




Figure 6. If Nonlinear Rotor/Propulsion System With Collective 

Input, Body Pitch Rate 
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Figure 6.1g Nonlinear Rotor/Propulsion System With Collective 

Input, Body Yaw Rate 
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Figure 6.1i Nonlinear Rotor/Propulsion System With Collective 

Input, Body Downward Velocity 




6.2 COMBINED PERTURBATIONAL AND GENHEL SIMULATIONS 


The perturbational rotorcraft/propulsion model was appended 
to the GENHEL simulation. The perturbational model consisted of 
the propulsion system constructed in Chapter IV and the 
rotor/fuselage system constructed in Chapter V. The 
perturbational model was connected so, at the option of the 
person running the program, the perturbational rotor, fuselage 
aerodynamics, and propulsion subroutines would replace those of 
the GENHEL simulation. The organization of the resulting model 
is shown in Figure 6.2. Note that there are two possible paths. 
The first is through the existing GENHEL nonlinear rotor and 
attached propulsion subroutines, and the second is through the 
perturbational subroutines. 

Recall that the perturbational equations do not yield total 
values. The perturbational equations are integrated separately 
from the GENHEL equations. After the perturbational equations 
are integrated, trim values are added to the perturbation values 
to yield total values which can be compared with the GENHEL 
predictions . 

Figure 6.3 compares the responses predicted by the 
perturbational equations with those of the GENHEL simulation. 

The trim condition here is the same one used previously; 120 
knots forward speed, 1000 feet altitude. The solid lines in the 
figures are the GENHEL results and the dashed lines are the 
perturbational model results. Note that the problems caused by a 
poor fit of the rotor equations (this was discussed in the 
previous chapter) appear here also. 

6.3 DISCUSSION 

A mostly linear perturbational model for modeling the 
interactions between the rotor system and the propulsion system 
of a helicopter was developed and demonstrated. The program that 
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Figure 6.2 Organization of Coupled GENHEL and Perturbational Simulation 
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implements this model has a modular structure that provides 
flexibility to change, add or drop independent modules from 
consideration on any one run. 

The perturbational model has been coupled with a nonlinear 
rotorcraft simulation. This coupling provides a means for 
checking out the model (other than with real test data) and 
provides a flight control system to "fly" the model. 

The input-output software and plotting software developed 
for the nonlinear simulation are then also available to the 
perturbational model. 

One rotor/propulsion interaction problem, rotor droop was 
demonstrated with the perturbational model. The model was also 
developed to look at other problems, such as hunting between 
engines, but this application was not tested. 

The generic model formulation methodology checked out, but 
the model (consisting of the coefficients) did not. This was 
because of the difficulty encountered in identifying a rotor 
model. The periodic effects, as discussed in Appendix A, were 
not sufficiently understood before creating the data file for 
identifying the rotor model and performing the OSR runs. Since 
then, more has been learned about the periodic effects, most 
importantly their structure in multi-blade coordinates. The 
rotor model that was obtained was sufficient to test and check 
out the generic model methodology, but the coefficients did not 
match the response of the nonlinear rotor model. 
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VII. SUMMARY AND RECOMMENDATIONS 


7.1 MODELS DEVELOPED 


SCT has developed several models during the period of this 
contract. These models include: 

(1) a nonlinear rotorcraft propulsion system model 
including two engine subroutines, two fuel control 
subroutines and a drive train subroutine. 

(2) a mostly linear propulsion system model including 
two linear engine models, two nonlinear fuel 
control models and a linear drive train model, 

(3) a mostly linear, perturbational rotor system model 
including a linear perturbational rotor model in 
multi-blade coordinates, a perturbational fuselage 
aerodynamics model in multi-blade coordinates, and 
a nonlinear lag damper model, 

(4) an augmented version of the GENHEL simulation 
which consists of GENHEL [14] augmented to include 
the nonlinear propulsion system in (1) above, 

(5) a perturbational rotor/propulsion model, which has 
been connected to the augmented GENHEL simulation 
in (4) so it can make use of the GENHEL flight 
control system, and which has been parameterized 
to model the RSRA. 

The nonlinear and the mostly linear propulsion system models 
were developed as modules that could be connected into a 
rotorcraft system. The perturbational rotor and fuselage 
aerodynamics models were developed similarly as modules. The 
perturbational rotorcraft model was constructed by connecting the 
mostly linear propulsion system model and the perturbational 
rotor and fuselage aerodynamics models together. 

Some of the models developed during this contract were 
checked out by comparing their responses with the responses of 
existing nonlinear simulations. No test data was analyzed. The 
drive train model, on the other hand, could not be checked out in 
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this manner because no previous drive train model was available 
with parameters for comparison. Further work is necessary to 
check out this model. 

The perturbational fuselage aerodynamics model and the 
perturbat ional rotor model outputs were compared to outputs from 
the nonlinear GENHEL simulation for verification. The fit for 
this perturbational rotor model has been poor. The poor fit was 
attributed to coefficients in the model that were in error 
because the periodic effects in the rotor dynamics were not 
considered when the coefficients were identified from GENHEL 
response data. Analyses of the periodic effects in this and an 
associated study [17] have since determined a better form and 
approach for modeling these periodic effects in multi-blade 
coordinates. Further work should be performed to identify new 
rotor model coefficients making use of this new information. 

7.2 RECOMMENDATIONS 

Recommendations for further work with the perturbational 
rotorcraft models and generic methodology developed during this 
contract include the following: 

• Develop parameters for a new perturbational rotor 
model. Include the effects of the periodic 
dynamics as described in Appendix A. 

• Test the perturbational model in control system 
design . 

• Analyze wind tunnel or flight test data by using 
the perturbational model as an equation structure 
for parameter identification and as a simulation 
for executing and testing the identified model. 
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APPENDIX A 

TRANSFORMATION FROM ROTATING TO MULTI -BLADE COORDINATES 


A.l MULTI -BLADE TRANSFORMATION 

The transformation from rotating to multi-blade coordinates 
IS defined, for flapping of a five bladed rotor system, by 
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multi-blade lagging coordinates 
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^r" ^2’ ^4’ ^5^ 

= rotating frame lagging angles 

T 

H = same as above 


The transformation, H, has the interesting properties 

H’^= I (A. 4) 

H * H J 



= main rotor speed 


Now, the equations of motion for the flapping and lagging of 
the blades are given in the rotating system by 




control inputs 
and fuselage 
inputs (A. 5) 


where, assuming the blades are independent and identical, the 
coefficient matrices, , . . . , ) , are 

diagonal. For example. 
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But, the multi-blade coordinates 
individual blade angles (6^. and 
system via the transformation 
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which follows from Eqs. (A.l) to (A. 4). Thus, applying this 
transformation (Eq. A. 6) to the blade equations in rotating 
coordinates (Eq. A. 5) yields the blade equations in multi-blade 
coordinates . 
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Figure A.l Example Dynamics Matrix Structure for Five-Bladed Rotor 
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for the generic model, the terms H^FH, H^F'H, etc. are 
diagonal and given by, for example; 

h'^FH = fQ [I] (A. 10) 

where [I] = S x 5 identity matrix 

When periodic effects are not small, then there are 
additional coefficients on and off the diagonal of the "constant" 
matrix, as well as matrices of coefficients for periodic terms 
(e.g. the cos 5 term in Figure A.l). 

A useful theorem about the solution of a periodically time 
varying equation at trim is found in [19]. This theorem states 
that the trim solution of a periodically time varying system at 
trim contains a constant part and a periodically time varying 
part . 

^trim “ ^const * ^(t) periodic (A. 11) 

The trim value is not constant but time varying. 

A. 2 IDENTIFICATION FORM 

The terms in the equations of motion (Eq. A. 8) for the rotor 
in multi-blade coordinates can be separated into several groups, 
as was done in Ref. 20. First, there is no need to include the 
portions of the blade accelerations containing the "J" terms in 
the OSR regression analysis, since they are already known. 

These "J" terms (i.e. 2J, j, J^) can be subtracted off 
before the regression analysis. Second, as was found in the 
present study, the diagonal form of Eq. (A. 10) can be assumed to 
simplify the regression for a constant coefficient model. In 
addition. Ref. 16 shows that if the diagonal structure of Eq. 

(5.10) is not imposed, the mode shapes cannot be correctly 
predicted by a constant coefficient model. Finally, the fuselage 
input portion of Eq. (A. 8) can be studied to determine "known" 
terms such as inertial effects which can also be subtracted from 
a regression analysis. 
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Rearranging and multiplying out yields the useful form 
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The terras such as h'^FH (also h"^GH, h'^G’H, h'^AH, 

H^A'H, etc.) have an interesting forra, as shown in Refs. 16 and 
17, and explained below. If the coefficient matrices ^ ) > 

^*('j^)» are expanded in series and 

truncated, for example, at second order, 

F(^) = Fq + Fj^costl' + F 2 Q cos24' + Fj^g siml^ 

+ p 2 g sin24' (A. 9) 

then, after much trigonometric algebra, we get the form shown in 
Figure A.l. 

The fo terms are the components that would appear in a 
constant coefficient model (i.e., one which neglects periodic 
rotor dynamics) [16,17]. The component ^ ^lc»^2c» 
associated with the periodic effects. For small advance ratios 
(u <.3) [18], the periodic effects, and thus these coefficients, 
are small. 

The importance of these results to the present study is 
that, for a constant coefficient model such as is being developed 
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(2) defined a methodology for parameterizing the model structure to 
represent a particular rotorcraft; (3) construction of a nonlinear coupled 
rotor /propuls ion model as a test case to use in analyzing coupled system 
dynamics; and (4) an attempt to develop a mostly linear coupled model 
derived from the complete nonlinear simulations. Volume I contains 
the details of the modelling process and its implementation approach. 
Volume II contains documentation of the computer models developed under 
this investigation. 
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